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We derive quantum kinetic equations for fermion and bo- 
son production starting from a (j> 4 Lagrangian with minimal 
coupling to fermions. Decomposing the scalar field into a 
mean-field part and fluctuations we obtain spontaneous pair 
creation driven by a self-interacting strong background field. 
The produced fermion and boson pairs are self-consistently 
coupled. Consequently back reactions arise from fermion 
and boson currents determining the time dependent self- 
interacting background mean-field. We explore the numerical 
solution in flux tube geometry for the time evolution of the 
mean-field as well as for the number- and energy densities 
for fermions and bosons. We find that after a characteristic 
time all energy is converted from the background mean-field 
to particle creation. Applying this general approach to the 
production of "quarks" and "gluons" a typical time scale for 
the collapse of the flux tube is 1.5 fm/c. 
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I. INTRODUCTION 

An ultra-relativistic heavy ion collision is a very complex 
process and its theoretical description requires a fundamen- 
tal understanding of all phases M . After the collision of the 
nuclei, it is expected that a strongly coupled quark gluon 
plasma is formed Q. The production of hard partons due 
to collisions is one successful way to describe properties of 
the pre-equilibrium phase Alternatively the flux-tube 

model was developed describing the spontaneous creation of 
soft particles in a chromo-electric background field Jj|-fic|. 

The idea is based on the QED picture of spontaneous vac- 
uum pair creation in a strong field ^],[l| ■ The strong electric 
background field leads to a restructuring of the Dirac sea and 
drives the system into a false, unstable vacuum which decays 
by emitting particles. Experimentally this was never mea- 
sured because of the insufficiently small fields compared to 
the mass of the electrons in e.g. optical laser experiments. 
Planned new facilities such as the X-ray free electron laser 
(XFEL) may be able to probe the domain of non-perturbative 
QED and melt the vacuum . 

In QCD the string tension is large compared to the mass of 
the partons. Therefore the main application of the Schwinger 
mechanism is for ultra-relativistic heavy ion collisions to study 
the production of a strongly coupled system of partons within 
the flux tube model Different versions of a Schwinger-like 
source term have been suggested, e.g. fll5j-|l7|] , and are imple- 
mented in transport equations to describe formation and evo- 
lution of a quark-gluon plasma on a phenomenological level. 



However the precise connection between field theoretical 
approaches and a kinetic theory is very challenging and an 
unsolved problem in general. For the considered example of 
the Dirac/Klein-Gordon equation in an external field this con- 
nection was found in The resulting quantum Vlasov 
equation was successfully applied to fermion- and boson pair 
creation 

In the applications considered so far either fermions or 
bosons were treated. One important aspect of most of the 
studies was the appearance of plasma oscillations due to back 
reactions ^i|j22| and their damping due to collisions po| , p3|l . 
Herein we introduce a self-consistent scheme in which fermion 
and boson fields are coupled via a strong background mean 
field. This internal field replaces the external field description 
of recent approaches and dynamically couples fermion and 
boson production. This coupling is most important for the 
back reactions since the total induced current consists of both 
fermionic and bosonic contributions. We solve the result- 
ing coupled equations in cylinder geometry j24| to implement 
physical boundary conditions leading to a lower bound in the 
infra-red momentum region. This procedure not only satis- 
fies the geometry constraints characteristic for a flux tube but 
also allows the production of massless scalar particles ("glu- 
ons" ) . The fermionic "quark" mass is treated as a momentum 
and temperature dependent quantity as obtained from QCD 
Dyson-Schwinger equation investigations |l7|]2l| . 

The article is organized as follows. In Section II we intro- 
duce the effective Lagrangian which we use to derive the ki- 
netic equations. In Section III we introduce cylindric bound- 
ary conditions for which we solve the equations numerically 
in Section IV. We summarize our results in Section V. 



II. EFFECTIVE LAGRANGIAN AND EQUATION 
OF MOTIONS 

In this article we will adopt the flux tube picture. One 
important aspect is the particular geometry assuming that the 
"chromo-electric" field acts along the longitudinal direction. 
We introduce the arbitrary fixed 4- vector n M which allows the 
implementation of the flux tube geometry as described below. 

We start our discussion from the following Lagrangian for 
fermions and bosons: 
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L{x) = d^*(x)8^{x) - m 2 + |$(:r)| 2 - e 2 |$(x)| 4 
+ie+ {[V„$*(:r)]$ 2 (:r) - [$*(a;)] 2 V n $(x)} 



($(a;) + <3>*(:r)) %l!(x)( / yn)ip(x) 



(1) 



where (-) denotes fermions, ip, and (+) bosons, <J>. The op- 
erator V n = permits to introduce the self-interacting 
scalar field with first order derivative only. The total field, 
$(a;), can be decomposed into a mean field contribution, $0, 
and fluctuations around it, ip, (e.g. JM) 



(2) 



We consider $o(i) as a neutral, space-homogeneous back- 
ground mean field: $0 =< $ >■ The field of the fluctuations 
tp{x) is in general complex corresponding to a charged field 
with a vanishing mean value < tp(x) >= 0. The Lagrangian 
(0) contains in particular terms of the order |<E>(:r)| 4 related 
to self-interaction. As we will discuss below, we keep the 
^(a;)! 4 contribution to determine only the background field 
<£{)(£)■ That makes possible a proper inclusion of back reac- 
tions. However we neglect all higher orders in the fluctuations 
in deriving the equations of motion and the kinetic equations. 
That means we restrict ourselves to the limit without colli- 
sions.^] We obtain the following Lagrangian for the quantum 
fluctuations in mean field approximation 

L ~( x ) = ^V'( a; )7 M <9 M i>( x ) - m-ip(x)ip(x) 

+e-<$> {t)^{x){n^)%}j(x) , (3) 
L + (x) = d^,ip* (x)d tl ip(x) — m\ \ip{x)\ 2 

+ie+$ (i) [(V„<p*(x))<p(x) - ip*(x)V„ip(x)] 

-4l<Mt)| 2 |^)| 2 • (4) 

The corresponding equations of motion are found within stan- 
dard techniques and read 



iy^d/i — e_7 M n M $°(f) — m_ \ip{x) = 



[D*^ + m 2 + ] <p(a 



0, 



fiS*o(t) + m 2 + Mt) + 44*g(t) - 3- -j+=0, 



(5) 

(6) 
(7) 



where Eq. (|5|) is the Dirac equation, Eq. (|s|) is the Klein- 
Gordon equation with the covariant derivative 



D u 



(8) 



and Eq. (^) provides the time dependent background mean 
field. It is important to observe that the external gauge field 
which typically appears in applications for pair creation in 
strong fields in the Dirac and Klein-Gordon equation is re- 
placed by an internal mean field $0. The third term on the 
left-hand side of Eq. (mj is non-linear and appears as result 



of the |$| 4 self-interacting contribution in the Lagrangian (|l|). 
The second (mass term) and the third term are not contained 
in the external field description]^] The equations of motion 
are self-consistently coupled, i.e. the mean background field 
acts on the fluctuations and vice versa: the produced charged 
particles generate the currents forming the background field. 
That is the well-know back reaction phenomenon. The con- 
traction of these currents with n M reads 

j-(t) = -e_ < ii(x)('yn)tp(x) >, (9) 
j+(t) = -e+ < if* (x)V „<p(x) - i(X7 n ip*(x))tp(x) 

-2e + $ (t)ip*(x)(p(x) > . (10) 

In the following we will identify the fermions as "quarks" and 
the bosons as "gluons", i.e. we consider the "gluons" to be 
scalar. Alternatively one could consider the bosons as a vector 
field and such an approach is demonstrated in the Appendix 
A yielding the same Lagrangian as Eq. (^). Many aspects of 
QCD are not contained in our approach and clearly to identify 
the fermionic and bosonic degrees of freedom with "quarks 
and gluons" is a very optimistic view. However we hope that 
some qualitative features discussed below are robust and also 
hold for more sophisticated theories than our toy model. 

One way to proceed is to solve the equations of motions 
directly. We are interested in finding kinetic equations being 
the exact analogue to these equations of motion. Starting 
from Eqs. (^) and (H) it is possible to derive such quantum 
kinetic equations. In |l^ , p^| it was explained in detail how to 
introduce quasi-particles, diagonalize the interaction Hamil- 
tonian and describe the transition from an unstable to a stable 
vacuum as a dynamical process employing a time-dependent 
Bogolyubov transformation. The resulting equation for the 
single particle distribution function is exact on the mean field 
level in a space-homogeneous field, i.e. it preserves the quan- 
tum statistical nature and the pair creation phenomenon. We 
choose n M = (0, 0, 0, 1) which is equivalent to considering a 
vector potential in temporal gauge with a background field 
acting in z-direction. We obtain 



df±(p,t) , M df±(p, t) 

L e±a(t) — = b±(p,t), 



dt 



dpu 



(11) 



where a(t) — —d^o(i)/dt is the space-homogeneous time de- 
pendent "chromo-electric" field strength, f±(p, t) is the distri- 
bution function of the partons in quasi-particle approximation 
and S± (p, t) is the source term 



S±(p,t) = ~W±(p,t) I dt'W±(p,t,t') x 



[l±2/±(P± (*,*'),*')] cos J 2 J druj ± (p;t,r) \ (12) 



1 A systematic way of defining the mean-field approximation 
and beyond is the large N expansion, see e.g. |27[|. 



2 Note that for a particular Lagrangian of similar struc- 
ture, Eq. ([?]) has a model specific shape, e.g. starting 
from a Ua(1) symmetry breaking Lagrangian of the Witten- 
DiVecchia-Veneziano type the non-linear contribution can be 
a transcendent function |^,^8] . 
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describing the creation and annihilation processes of fermions 
and bosons as result of vacuum tunneling in a strong quasi- 
classical mean- field $o(£)- The transition amplitudes are 



W±(p;t,t') 



e ± a{t')P±(t,t') 
<4(p;M') 



£±J 



P±(M') 



2s± 



(13) 



For fermions s_ = 1/2 and for bosons s+ = 0. The energy 
squared of the quasi-particles reads 



u±(p,t,t') =£±±+P±(t,t'), 
t 

P± (*,*') =p\\ ~ e ±y dra(r) 



(14) 
(15) 



where e^_ x = m± + Pi + p| is the transverse energy. It is 
assumed that lim a(t) = and we simplify the notation by 

t — * — oo 

introducing 



w±(p) = w±(p,t,t) 



W ± (p,t) = W±(p,t,t). 



(16) 



The kinetic equations (^lj) coincide with recent results and 
retain all physical information contained in the field equa- 
tions (|^)-@- The main properties of the source term are 
the inclusion of quantum statistical effects due to the Fermi 
suppression (Bose enhancement) factor for fermions (bosons) 
and the nonlocal time structure. "Quarks and gluons" are 
produced with a non-trivial, non-equilibrium momentum dis- 
tribution which is different in transverse- and in longitudinal 
direction to which the field $o couples according to the cho- 
sen anisotropy. Apart from the conceptual interest of a ki- 
netic formulation, the practical merits are obvious: e.g. the 
straightforward identification of the Markovian- or low den- 
sity limit and the phenomenologically simple inclusion of col- 
lisions. 

The background field $o decays through particle creation 
and generation of currents which are readily obtained: 



j±(t) = 2e±g± 



d 3 p 



(2 7r) 3 \u±(p) 
w±(p) df±(p,t) e±&(t)p\\ 



f±(P,t) + 



e±a(t) dt 



<LL± 

P\\ 



2s-)- 



(17) 



where g± is the degeneracy factor The first term is propor- 
tional to the occupation number itself and the second term 
is related to polarization being proportional to the rate. The 
last term is a regularizing counter term. Its origin is associ- 
ated with charge regularization; further details can be found 
in @|. 



III. CYLINDRIC BOUNDARY AND THE FLUX 
TUBE GEOMETRY 



finite volume. We introduce this feature because of two rea- 
sons: (i) to implement the given anisotropic symmetry and 
(ii) to avoid difficulties with massless particle production for 
"gluons". It is clear that a finite size in transverse direction 
leads to a quantization of the transverse momentum, so that 
at zero mass m+ = the minimal energy at which particles 
can be produced is A IR . 

One possible, simple way is to introduce the boundary con- 
ditions 



<&{X3,P0,<t>,t) = 



(18) 



to require that the field vanishes outside a fixed flux tube of 
the radius po, e.g. [Q. The Klein-Gordon equation in such 
cylinder coordinates reads 



dl - vi - dl 



(p(x 3 ,p,(p,t) = 0, 



(19) 



Id Id 2 
+ 



(20) 



where the Laplace operator is given by 

Vi= * 

dp 2 p dp p^ 

The radius is denoted with p and the angle is <f>. Here we 
assume that the mean field depends only on time and not 
explicitly on the radius: $o = < I ) o(i)J^ This assumption allows 
to make a separable ansatz of the following form 

<p{x 3i p,4>,t) = Ne a +e~ ih >**T(t)R(p). (21) 

Substituting ( |2l| ) into ( |l9| ) we find two equations, one for the 
time dependent part 

f(t) + [m 2 + + (jka - e+$o(t)) 2 + e 2 } T(t) = (22) 

and one for the spatial dependence 

P 2 R (p) + P R (p) + (e 2 p 2 - I 2 ) R(p) = , (23) 

where e is any constant value. The latter equation is known 
as the Bessel equation and with the boundary condition 
R(epo) ~ 0. Eq. ( ^lj ) can be rewritten in terms of Bessel 
functions 



(Pnl,k 3 (x3,p,<t>,t) 



-ik a x arr m Jl{£nlp) 



v^poJ; {e n ipo) 
where e n ipo is the n — th zero of the Bessel function J;. 



(24) 



po 


1 fm 


2 fm 


3 fm 


A 1H 


0.48 GeV 


0.24 GeV 


0.16 GeV 



TABLE I. The radius of the flux tube po is fixed and leads 
to a discretized momentum with a lower bound A IR . 



Another element that we would like to introduce in our 
model is a cylindrically confined region for the particle pro- 
duction process. A similar approach was explored in [MM] 
and an advanced description of a dynamical flux tube was 
given in Q3o|| . The background field is constructed to act only 
inside a cylindric tube and is assumed to vanish outside this 



' j A p- dependent field $0 would lead to nontrivial changes 
in the quasi-particle representation itself and to more compli- 
cated final expressions for the equations of motion. 
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The resulting kinetic equations are identical to Eqs. 
but the quasi-particle energy is now discrete and reads 



lip, t) 



+ (fc 3 -e±$ (i)) 



2 2 2 

£±_Lnl = m ± + £nl ■ 



(25) 
(26) 



The kinetic momentum pu — kz — e±&o(t) contains the time 
dependent mean-field and the canonical momentum kz. The 
lowest possible values for the transverse energy are assumed 
for Eio and the corresponding A IR — eio are given in Table 
I for three different flux tube radii. This means that the 
transverse energy never vanishes. Even for a zero gluon mass 
the discretization provides a lower bound for the momentum. 
To proceed we consider all observable values as mean values 
in the following way 



<< 



>>P= 2 / P d P 

n Po Jo JO 



< ... > 



(27) 



Assuming such mean (p- independent) quantities is a very 
simple approximation. However for the qualitative study per- 
formed it is sufficient and can certainly be improved for a 
particular experimental application. The charges in strong 
coupling limit are fixed to be e± = 1 throughout the numeri- 
cal calculations and the degeneracy factors for a three flavor 
system are for quarks g- — 18 and for gluons g+ — 16. Quan- 
titatively the numerical results depend on these values, qual- 
itatively the main effects discussed below are robust against 
small changes of all parameter values. 




0,25 0,5 
p [GeV] 

FIG. 1. Momentum dependence of the quark mass for 
different temperatures obtained from solving the QCD 
Dyson-Schwinger equation for the quarks employing a sim- 
plified gluon propagator. The fixed radius of the flux tube, 



p = 2 fin, 
GeV. 



leads to a lower momentum bound of A — 0.24 



Before we present our numerical results we briefly dis- 
cuss the mass scales of the produced particles. The "gluon" 
mass is chosen to be zero or finite. This is possible because 
of the geometry constraint discussed above. The fermion 
mass is very small on the hadronic scale in the perturbative 
regime, m_(p — > oo) << Aqcd- In the non-perturbative 
domain dynamical chiral symmetry breaking appears and 
the quark mass increases by about two orders of magnitude 



m-(p — > 0) ~ Aqcd- For very large temperatures the quarks 
can be well approximated by their current quark mass for 
all momenta. However in vicinity of the phase transition a 
dynamical mass is generated. Furthermore, non-perturbative 
effects are manifest up to 2 — 3 x T c ~ 300 — 450 MeV. It was 
shown in a variety of model approaches that within the parton 
creation model employed in our investigation such tempera- 
tures are initially reached. Due to a rapid expansion the tem- 
perature will further decrease and non-perturbative dressing 
of the quarks appears. 

In QCD the momentum dependence of the quark mass can 
be calculated on the lattice, e.g. in |3lf i n Landau gauge. 
QCD Dyson-Schwinger equation models [ [l7p^ , p2] ] are in good 
agreement with these studies and we employ a simple model 
introduced in j3^| extended to non-zero temperature and den- 
sity in I17U34I . The instantaneous version of such an infrared 
dominant model Jn],^ leads to the following QCD gap equa- 
tion for the scalar part of the self-energy: 



b(p) = m + r\ 



m-{p) 



y/p 2 + m 2 _(p) 
[l-f e J(p,T,n)-f°J(p,T,p)] 



(28) 



and the vector part 



a(p) 



26(p) 
m + b(p) 



(29) 



where m_ = b/a. fZ q and /f 9 are the Fermi distribution 
functions for particles and antiparticles, respectively. The 
value of the mass scale parameter 7] = 1 GeV is inspired by 
the potential energy in a QCD string at the confinement dis- 
tance, V q q(r = lfm) ~ 1 GeV. The momentum dependence 
calculated with Eq. ( pj| ) is depicted in Fig. 1 for different 
temperatures. The chemical potential is fi — and the cur- 
rent quark mass mo = 5 MeV. For further calculations we will 
use the shape corresponding to T = T c ~ 170 MeV. For the 
current application a constant fermion mass with the value of 
mo would provide a good guide. However for future applica- 
tions, it will be very useful to have its complete temperature 
and momentum dependence introduced already. 



IV. NUMERICAL RESULTS 

The final equations solved numerically are given in the Ap- 
pendix B. In the following we compare the model cases (i) 
coupled fermion and boson production with non-linear self- 
interaction included; (ii) only boson production with self- 
interaction, i.e neglecting the fermionic contribution, j_ in 
Eq. (Q) to (iii) boson production without self-interaction, i.e. 
neglecting the nonlinear term in Eq. (Q) and j- . Further- 
more we will explore the dependence of the numerical results 
on the parameters of the model, i.e. the masses and the flux 
tube radius, see Table I. We use (a) a zero "gluonic" mass 
?n+ = and compare the results with (b) a finite "gluon" 
mass of m+ = 0.5 GeV; the mass of the quarks is given by 
the simple model defined by Eq. (Eq). 
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A. The mean background field 

We start the discussion of the numerical results with the 
background mean field er(t). As initial condition for the field 
we choose e+a(t = 0) = 4 GeV 2 corresponding to a large 
initial energy density of about e(t = 0) ~ 700 GeV/fm 3 re- 
lated to future LHC experiments. The time dependence of 
the "chromo-electric" field strength a(t) is plotted Fig. ti. 



Gluons 
Quarks 

Gluons and Quarks 



> 
O 



to 

+ 




t [fm/c] 

FIG. 2. The time dependence of the strong background 
"chromo-electric" field, cr(t), for "quarks and gluons" with 
self-interaction. The flux tube radius is po = 2 fm and the so- 
lutions for two different "gluon" masses m+ = and m+ =0.5 
GeV are compared . For zero "gluon" mass the flux tube col- 
lapses at Tf ~ 1.5 — 2 fm/c and all energy is transformed from 
the mean field to the created particles. For massive "gluons" 
weak plasma oscillations survive. 

The self-interaction leads to oscillations on a time scale of 
T p i ~ 0.5 fm/c in the early phase of the evolution. This re- 
sult is basically independent of the masses used, however it 
strongly depends on the initial field strength. The amplitude 
of the oscillations is damped due to back reactions and the 
background mean-field vanishes after a few periods at about 
Tf ~ 3r p ; ~ 1.5 fm/c for the case of zero gluonic mass. With 
other words: the field leading to pair creation disappears, the 
flux tube dynamically collapses. For m+ = 0.5 GeV the oscil- 
lations at the beginning of the evolution are damped as well 
but at Tf the field does not disappear completely. It rather 
evolves into a permanently oscillating field with an amplitude 
of 10% of the initial value.^j This behavior for t > Tf can 
be compared to typical results in an external field approach 
where a constant background field or an impulse shape was 
applied and undamped plasma oscillations appeared. The 
production of quarks is much less effected, see Fig ^| and a 
physical explanation of this effect seems to be straightforward: 



> 
u 
O 








4 We suppose that a detailed study of this effect in the ap- 
plication presented in plpq] could lead to a similar result. 



t [fm/c] 

FIG. 3. The time dependence of the strong background 
"chromo-electric" field, cr(t), for a flux tube radius of po = 2 
fm for "quarks and gluons" with self-interaction (solid line), 
"gluons" only with self-interaction (dashed line) and "quarks" 
only (dotted line). The "quarks" are much less effected by 
the dynamical coupling compared to the "gluons", cf. Fig. 
^| Inclusion of the "quarks" leads to a faster damping of the 
oscillations in the beginning of the evolution. 

The field energy is transformed into the energy of the pro- 
duced partons, i.e. the amplitude of the oscillations for t < t/ 
is decreasing with time. In the vicinity of large magnitudes of 
<j(t) the current is zero since the partons stop their collective 
motion within the flux tube and at this time the distribu- 
tion function is symmetrically peaked around zero momen- 
tum For fermions Pauli blocking is most efficient and 
no particles can be produced in already occupied states. For 
bosons, however, the enhancement of particle production has 
its maximum at this time and a large number of bosons is 
produced. Therefore a direct flow of energy from the field to 
the gluons occurs, i.e. the damping of the mean field happens 
very sufficiently. With other words, the quantum statistical 
nature leads to a fast energy transformation for bosons but 
to a suppressed one for fermions. 

The origin of the collapse of the flux tube is the back reac- 
tion phenomenon |3(|, the time scale on which this appears 
is given by self-interaction, see Fig. ^. This becomes clear 
by neglecting the non-linear terms in Eq. (Q). That is one 
new result of the present investigation. The plasma oscilla- 
tions are smoothly damped only, i.e. the energy is slowly 
converted from the mean field to the produced particles. In 
that case the background vanishes at much larger times. 

There is only a weak dependence of the results for the 
mean-field on the flux tube radius. Using the different radii 
given in Table I, the time scale 77 is changed by about 10 
% but all qualitative effects are unchanged. For values of 
A IR ~ m+, the lower bound in momentum acts like a mass 
scale itself and leads to similar results as discussed above in 
connection with the finite gluon mass. Numerical examples 
of that consideration are given in the next subsection. 
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FIG. 4. The time dependence of the strong background 
"chromo-electric" field, a(t), for "gluons" with (solid line) and 
without (dotted line) self-interaction. The flux tube radius is 
po = 2 fm and m+ = 0. It is apparent that the inclusion 
of self-interaction leads to a rapid collapse of the mean back- 
ground field. Without self-interaction the energy is smoothly 
converted from the field to the particles. 



B. Particle production 

The self-consistent solution of the quantum kinetic equa- 
tions for quarks and gluons provides particle creation and a 
strong background mean-field discussed in the previous sec- 
tion. The large strength of the field initiates the particle 
production process and "quark"- and "gluon" pairs are spon- 
taneously produced. In Fig. [B| we plot the number density 



(30) 



of the produced "quarks and gluons". It is simple to verify 
that Eq. (^) is a finite expression, e.g. by using the equations 
given in Appendix B. At the beginning of the evolution the 
background field is very large and therefore many particles are 
produced in a short time. At the same time the background 
field decreases, energy is converted from the field into the cre- 
ated particles. At r/ ~ 1.5 fm/c the background mean-field 
vanishes and consequently no more quarks and gluons are pro- 
duced and the number density assumes a constant value. The 
number density is very large and therefore it is necessary to 
solve the full non-Markovian equation (|ll|). The application 
of the low density limit would underestimate the production 
of "gluons" and overestimate the "quark" production. Note 
that, in a realistic description of the further evolution of the 
produced plasma, collisions between the partons could equi- 
librate the system. Additionally, the number density would 
decrease when an expansion is implemented. 
The regularized energy density is defined by 



e±(t) 



g± 



E 



2tt 



w±(# (/±(P,*)-/±(iM)) (31) 
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c 
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Gluons p = 1 fm 
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FIG. 5. The time dependence of the number density for 
"quarks and gluons". A fast increase of the particle number 
is observed at the beginning of the evolution. At 77 ~ 1.5 
fm/c the background mean-field vanishes and therefore the 
creation process stops and the number density is constant. 
Different flux tube radii do not alter the results significantly. 



with the counter term |16 



f±(P,t) = 



e±a(t)p\\ 
4c4(p) 



P\\ 



2s±-i 2 



(32) 



The time dependence of the energy density is plotted in Fig. 
^| Due to the strong background mean-field many particles 
are produced at small times and energy is converted from the 
field to the particles (see Fig. ^). During the fast increase of 
the energy density a shoulder-like structure appears. The pe- 
riod of these waves is in tune with the oscillations of the mean 
field. The energy density reached at large times for "gluons" is 
larger compared to the "quarks" due to Pauli blocking acting 
on the quarks. The value of about 200 GeV/fm 3 corresponds 
to a (quasi-equilibrium) temperature of about 1 GeV and is 
characteristic for planned ultra-relativistic heavy-ion collision 
experiments. The insertion in Fig. [] shows the sensitivity to 
the "gluon" mass. A finite gluon mass leads to weak plasma 
oscillations in the background mean-field even for large times 
as discussed in connection with Fig. and therefore the en- 
ergy density oscillates in concert with these repeated creation 
and annihilation processes. 

Naturally the question arises whether the system is in equi- 
librium after the production of particles has stopped and no 
fields are imposed. Although the number- and energy density 
is constant for large times, the system is still in a strong non- 
equilibrium state. Further interactions such as inter-parton 
collisions could provide thermalization, but that is a challeng- 
ing theoretical question on its own. Further understanding 
can be provided by studying the pressure for the "quarks and 
gluons" given by the following expressions: 
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FIG. 6. The time dependence of the energy density for 
"quarks and gluons" for po = 2 fm. Due to the collapse of 
the background mean field, all energy is converted at 77 ~ 1.5 
fm/c into the partons and a constant value of the energy den- 
sities is assumed. Different "gluon" masses do not alter the 
results significantly, see insertion. 
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(35) 



(36) 



where the regularizing counter term f c is given in Eq. 
and v c reads 



v± (p, t) 



&4(g) 



Pll 



2sj- 



(37) 



In Fig. Q we plot the different components of the pressure. It 
is apparent that the considered system is still out-off equi- 
librium. It is interesting to observe that the difference is 
much more pronounced for "gluons" (lower panel) compared 
to "quarks" (upper panel) . The "gluon" production is unhin- 
dered and preferable in longitudinal direction, for "quarks" 
again Pauli blocking prevents a drastic difference. In an equi- 
librated system the longitudinal and the parallel pressure con- 
tributions would be equal. 




2 

t [fm/c] 

FIG. 7. The time dependence of the pressure of "gluons" 
(lower panel) and "quarks" (upper panel). The difference 
between the longitudinal and transverse part indicates that 
the system is out-off equilibrium also at large times. Pauli 
blocking prevents that the difference is less pronounced for 
"quarks" compared to "gluons". 



V. SUMMARY 

We have described coupled fermion and boson production 
within a quantum kinetic approach in cylinder geometry. The 
strong background field is given by a dynamical mean-field 
containing self-interactions. Back reactions are included, both 
fermionic and bosonic currents modify the initial "chromo- 
electric" background field. Strong particle creation appears 
at the beginning of the time evolution, however due to back 
reactions the background mean-field vanishes. The time scale 
at which such a collapse of the flux tube happens depends 
on the self-interaction. This general approach has been ap- 
plied to the production of " quarks" and "gluons" and using a 
zero "gluon" mass, a momentum dependent quark mass from 
Dyson-Schwinger studies in concert with a flux tube radius of 
about 2 fm, we find a typical time scale of 1.5 fm/c for this 
effect. At that time all energy is converted from the field to 
the particles and the evaluated number- and energy densities 
reach a constant value being of the typical order of magnitude 
of a few hundred GeV/fm 3 . The system is still out off equi- 
librium, we have exemplified this by calculating the pressure 
components which are different in longitudinal and transverse 
direction. 

The further evolution of the system will strongly depend 
on additionally included interactions, such as collisions. A 
detailed study will be reported elsewhere. 



ACKNOWLEDGMENT 

We thank CD. Roberts for helpful discussions. This 
work was supported by Deutsche Forschungsgemeinschaft 
under project number SCHM 1342/3-1, AL 279/3-3 and 
436 RUS 17/102/00; and partly by the Russian Federation's 



7 



State Committee for Higher Education under grant N E00- 
33-20. 



APPENDIX A: VECTOR INTERACTION AND 
BOSONS 

In Section II we have introduced an effective Lagrangian in 
mean field approximation based on coupled boson and fermion 
fields. The boson fields have been considered as scalar fields. 
In this appendix we will briefly demonstrate that starting 
from a charged vector field, the same Lagrangian in mean 
field approximation can be obtained. 

The fermionic contribution is identical to Eq. (|^) and there- 
fore restrict ourselves to the bosonic contribution. The La- 
grangian for a minimal substituted charged vector field reads: 



(Al) 



where D v — d v + ie g uo v is the covariant derivative. The de- 
composition into a mean-field contribution and fluctuations 
reads 



<j3h = &n + Sf,., 



(A2) 



Using this Ansatz in Eq. (Al) we obtain 



(D^(x)y (D»^(x)) = 

(0„n*0"fi" + d v 5ld v 8' i - e\ jfij 4 - e+ |n| : 

+ie+ [d v 5* t 5 v ^ - St^d'S 1 *] 

+ie+ [duSi*s v s> > - sts^n"] , 



(A3) 



where we have neglected odd orders in the fluctuations 8. 
Furthermore we focus on the collisionless limit, i.e. we do not 



consider e 4 



terms. 



Assuming that the fluctuations do not depend on the mean- 
field part and imposing the initial spatial anisotropy of the 
flux tube, we can write in Hartree approximation 



(5* k 8 k ) =3(M 2 ) , (s;s k ) = 8 lk (M 2 ) 



(A4) 



where i, k = 1, 2, 3. ip is the effective fluctuating scalar field. 
The mean-field reads 



tt M = (0,0,0,$o(t)) 



(A5) 



Employing this Hartree-Fock like approximation in temporal 
gauge we can construct the following modified Lagrangian 

L v = ?>dv<£i*d v <j> — 3m+|<p| 2 — 3e+ \(p\ 2 $o(i) 

+ie+<E>o(£) [d3<p*<£ - <fi*d3<p] ■ (A6) 

Eq. ( |A€^ can be generalized, using 83 — > V n , to obtain the 
Lagrangian Eq. (ft). 



APPENDIX B: EQUATIONS FOR THE 
NUMERICAL SOLUTION 



Eq. ( jl ip is an integro-differential equation. It can be re- 
expressed by introducing 



dt 



'W±(p,t,t')(l±2f±(P±(t, «'),*')) x 
siu2 I druj±(t,T)], (Bl) 
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cos[2 / dTLu±(t,r)], (B2) 



dt 



u±(p,t) 



v±(p,t) 



with the initial conditions u(p, 0) = v(p, 0) = 0, in which case 
we find 

dV±iP ' t] +e ± a(t) dV±iP > t) 



-W±(p,t)(l±2f±(p,t)} 

- 2w±(p)u±(£*), (B4) 

+ e ± <j{t) du t {P,t) = 2u ± (p)v ± (p,t). (B5) 
°P\\ 



dt 

du±(p,t) 
dt 



where the total energy is defined in Eq. ( |25[ ) for bosons and 
in Eq. (jlij) for fermions. The transition amplitudes are given 
by 



W±(p,t) = 



e±a{t)p« ( e± ±n i 



w±(p) \ P\\ 



2s-j- 



(B6) 



For the currents we obtain the following final expression 

dpu p« 



3±(t) = 2e ± g ± ^ V / 



2tt w±(p) 



f±(p,t) 



(B7) 
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